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ABSTRACT 


Pressure-sensitive  paint  measurement  on  a  transonic  compressor  rotor  required  the 
prior  development  of  phase-locked  cumulative  imaging  using  a  disk-rotor  driven  by  a  high- 
speed Hamilton-Standard  turbine  as  a  developmental  test-bed.   The  turbine  was  installed  in  a 
protective  housing  in  the  Gas  Dynamics  Laboratory  and  connected  to  an  8000  cu  ft,  300  psi 
air  supply.    An  hydraulic  pump  provided  bearing  lubrication.    A  once  per  revolution  trigger 
signal  was  produced  from  a  light-emitting  diode  and  PIN  photodiode  pair.   The  imaging 
system  consisted  of  an  intensified  CCD  video  camera  externally  triggered  by  the  1/rev  signal 
via  a  waveform  shaping  circuit  designed  for  the  present  application.   Images  were  captured  at 
camera  gate  speeds  calculated  to  eliminate  pixel  blur  and  image  integration  times  were  varied 
to  optimize  image  intensity  and  spatial  resolution.    Structural  and  modal  analyses  of  the  disk- 
rotor  were  conducted  and  a  simplified  numerical  model  of  the  flow  was  computed.   Ratioed, 
colored  images  were  produced  for  wheel  speeds  to  20,000  RPM.    The  effect  of  the  radially 
varying  stagnation  temperature  was  evident,  underscoring  the  importance  of  quantifying  and 
accounting  for  the  PSP  temperature  sensitivity  so  that  quantitative  pressure  data  may  be 
obtained.    Recommendations  for  a  follow-on  program  are  reported. 
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I.     INTRODUCTION 


With  increasing  emphasis  in  the  design  of  turbomachinery  on  three-dimensional 
numerical  computation,  validation  of  new  design  and  analysis  techniques  has  become  critical. 
Toward  that  goal,  the  development  of  a  pressure-sensitive  paint  (PSP)  measurement  technique 
has  been  initiated  at  the  Naval  Postgraduate  School  (NPS).  The  technique  is  to  be  applied  as 
a  new  diagnostic  in  an  experimental  appraisal  of  the  numerical  methods  used  in  the  design  of 
a  rotor  for  the  NPS  Turbopropulsion  Laboratory  transonic  compressor.    Following  a 
thorough  evaluation  of  the  numerical  methods  used  to  arrive  at  the  present  geometry, 
including  the  validation  of  computational  off-design  solutions,  new  designs  can  be  attempted 
which  optimize  benefits  from  blade  sweep  and  tip  clearance  flows.    Thus,  the  full 
experimental  validation  of  the  new  rotor  design  is  key  to  attaining  higher  performance 
designs  at  lower  test  and  evaluation  costs.   The  present  study,  building  on  experience  gained 
in  preceding  work,  involved  the  construction  and  use  of  a  turbine-driven  PSP  test-bed  to 
develop  the  measurement  technique.   PSP  images  were  successfully  acquired  from  a  rotor 
operating  at  sonic  tip  speeds,  and  problems  which  must  be  solved  before  measurements  can 
be  made  on  the  compressor  rotor  were  clearly  defined.    A  summary  of  previous  related 
research  at  NPS  and  a  brief  overview  of  the  principles  of  the  PSP  measurement  technique  are 
provided  as  an  introduction  to  the  present  study. 

A.        SUMMARY  OF  RESEARCH 

The  design  of  a  stage  for  the  NPS  transonic  compressor  test  rig  by  Sanger  (1996) 
used  a  quasi-three-dimensional  Euler  code  using  distributed  body  forces  to  model  viscous 
effects.   Final  stator  and  rotor  blade  aerodynamic  geometries  were  derived  iteratively  to  attain 
the  design  goals.     The  stage  was  machined  by  Turbocam  Inc.  to  coordinates  calculated  for 
the  unloaded  condition.   In  preparation  for  test  of  the  Sanger  stage,  a  study  by  Reid  (1995) 
evaluated  the  process  described  by  Sanger  and  returned  the  NPS  test  rig  to  operational 
condition. 

Concurrently,  the  commitment  to  PSP  measurement  at  NPS  was  initiated  by 
Seiv wright  (1996)  in  a  shock-boundary  layer  interaction  study  of  a  flat  plate  in  an  under- 
expanded  free-jet.   PSP  was  used  to  map  the  surface  pressure  field.   The  flat  plate  study 
enabled  the  application  of  PSP  in  a  rotational  experiment,  conducted  by  Varner  (1995). 
Elsewhere  Burns  et  al  (1995)  and  Hubner  et  al  (1996)  reported  using  PSP  in  the  rotating 
frame.   In  the  Varner  study,  a  flat-plate  disk-rotor  powered  by  an  electric  motor  was  used  to 
attempt  PSP  pressure  measurements  at  a  transonic  rotor  condition.    Imaging  requirements 


were  identified  and  a  phase-locked  circuit  was  developed  to  acquire  data.   Qualitative 
measurements  were  obtained  at  subsonic  rotor-tip  speeds  (13,000  RPM)  and  calibration 
studies  were  conducted,  providing  the  impetus  for  a  careful  evaluation  of  the  technique  and 
application  in  a  higher  speed  regime. 

B  .        THEORETICAL  BACKGROUND 

The  fundamental  mechanism  in  PSP  is  the  photoluminescent  characteristic  of  an 
activated  'probe'  molecule;  platinum  octaethylporphyrin  (PtEOP)  in  the  present  study.  The 
probe  molecule  achieves  activated  energy  states,  both  in  singlet  and  triplet  form,  through  the 
absorption  of  excitation  energy  from  ultraviolet  light  stimulation.   Photoluminescence  is  the 
radiant  energy  emitted  when  a  molecule  returns  from  an  excited  triplet  energy  state  to  a 
ground,  singlet  state.   This  process  is  quenched  when  collisions  occur  between  two  triplet-state 
molecules  with  opposite  electron  spin  orientation.     The  collision  promotes  a  non-radiating 
transition  of  the  probe  molecule  from  its  triplet  to  its  ground  state.     Therefore,  the 
photoluminescent  transition  and  molecular  collision  transition  processes  compete. 

PtEOP  is  predisposed  to  excite  to  a  triplet  state  when  activated,  therefore  tending  to 
phosphoresce  when  it  returns  to  its  ground  state.  However,  oxygen  in  the  atmosphere  occurs 
naturally  in  an  excited  triplet  state  with  an  opposite  electron  spin  orientation,  and  quenches 
the  phosphorescent  process  when  collisions  between  it  and  the  probe  molecule  occur.   The 
rate  at  which  these  collisions  occur  is  dependent  on  the  oxygen  partial  pressure,  which,  in 
turn,  is  dependent  on  atmospheric  pressure.    Stern- Volmer  related  the  rate  constants  of  the 
processes  involved  to  the  intensity  of  the  photoluminescence,  and  derived  an  inverse  linear 
relation  between  intensity  (I)  and  pressure  (P)  when  referenced  to  a  control.   Equation  1  gives 
this  relation,  showing  a  dependency  of  the  rate  constants  on  temperature  (T),  where  the 
subscript  zero  is  a  reference  (usually  ""wind-off")  condition. 

i  =  A(T)  +  B(T)|-  (*) 

o 

By  combining  the  PtEOP  in  a  oxygen-permeable  silicon  resin  binder,  application  to 
an  aerodynamic  surface  in  the  form  of  paint  allows  the  static  pressure  on  the  coated  surface  to 
be  determined  by  relating  the  intensity  of  the  reference  or  "wind-off  condition  to  the  flow 
or  "wind-on"  condition.    The  surface  temperature  must  be  constant  for  both  the  "wind-off 
and  "wind-on"  images  to  be  related  through  the  A  and  B  rate  constants  in  this  form. 
Seivwright  (1995)  and  Kavandi  (1990)  provide  thorough  discussion  of  the  mechanics  of 
phosphorescence  and  oxygen  quenching. 


C .        CONCEPT  AND  APPROACH 

In  the  present  study,  the  development  of  a  high-speed  test-bed  for  rotary  PSP 
measurements  and  refinement  of  the  image  acquisition  technique  enabled  PSP  results  to  be 
obtained  at  four  wheel  speeds.   The  installation  of  a  high-speed  air-driven  turbine  to  rotate  a 
dual-blade  disk-rotor  provided  the  ability  to  obtain  PSP  images  at  near  transonic  conditions. 
Mechanical  and  vibrational  characteristics  of  the  disk-rotor  were  addressed.   The  use  of  a 
phase-locked,  monostable-pulse  circuit  triggered  by  an  amplified  once-per-revolution  pulse, 
allowed  image  acquisition  with  a  charge-coupled-device  (CCD)  camera.   Pixel  blur  and  image 
integration  time  issues  were  solved  and  a  process  for  the  acquisition  of  images  at  optimum 
resolution  was  determined.    Circuitry  requirements  for  higher  speed  imaging  and  a  disk-rotor 
redesign  to  enable  higher  speeds  were  identified.   Most  importantly,  the  association  of  image 
intensities  to  radially  increasing  stagnation  temperature  was  verified. 

In  the  present  document,  Section  II  discusses  the  flowfield  produced  by  a  rotating 
disk  in  a  stationary  flow  and  presents  the  results  of  a  Navier-Stokes  numerical  computation. 
Section  III  details  the  experimental  setup  and  test  considerations.   Section  IV  outlines  the 
procedure  used  in  obtaining  four  high-speed,  ratioed  images  at  progressively  increased  rim 
speeds.   The  results  are  discussed  in   Section  V  and  conclusions  and  follow-on  research 
recommendations  are  given  in  Section  VI.     Appendices  A  through  E  provide  details  of  the 
numerical  code  used,  engineering  drawings,  postprocessing  software  scripting,  circuit 
specifications,  and  proposed  circuit  and  disk-rotor  redesigns. 


II.    FLOWFIELD  CHARACTERIZATION 

A .        OVERVIEW 

The  flow  created  by  a  finite,  high-speed  rotating  disk  within  an  otherwise  stationary 
environment  is  viscous,  compressible,  and  three-dimensional  (3D).   The  surface  friction 
forces  induce  radial  and  tangential  velocity  components  on  the  fluid  in  the  vicinity  of  the  disk 
surface.   For  continuity,  the  radially  outflowing  fluid  near  the  disk  is  fed  by  an  axially 
inflowing  fluid  mass  along  the  rotational  axis.   The  disk  is  therefore  a  centrifugal  type  pump. 
An  analytic  solution  for  the  laminar  boundary  layer  case  derived  by  Schlichting  (1960) 
assumed  incompressible  flow  and  zero  radial  pressure  gradient.    Calculating  Reynold's 
number  using  disk  radius  as  the  length  scale, 

Re  =  °*L  (2) 

where  co  is  rotational  frequency  and  "u  is  kinematic  viscosity,  Schlichting  showed  close 
agreement  between  analytic  and  experimental  displacement  thicknesses  for  both  laminar  and 
turbulent  conditions.  The  transition  from  laminar  to  turbulent  flow  was  shown  to  occur  at  a 
Reynold's  Number  of  approximately  3xl05. 

The  higher  disk  rotational  speeds  in  the  present  application  give  rise  to  a  turbulent 
boundary  layer.     To   account  for  compressibility  and  turbulence  and  not  be  constrained  to 
thin  viscous  layers,  a  numerical  solution  was  undertaken  using  the  external-flow  3D  Navier- 
Stokes  code  OVERFLOW.     The  computational  development  and  solutions  for  varying  disk 
rotational  frequencies  are  discussed  below. 


B  .        COMPUTATIONAL  DEVELOPMENT 

1 .        Code  Overview 

OVERFLOW  solves  the  three-dimensional  Navier-Stokes  equations  implicitly  using 
finite  differencing  and  block  tri-diagonal  solvers.  The  OVERFLOW  code  is  written  in  Fortran 
programming  language  for  workstation  use.  The  code  makes  use  of  an  input  file  which 
defines  the  time  step,  a  minimum  Courant-Freidrichs-Levy  (CFL)  number,  external  flow 
boundary  conditions  and  the  solution  variables.    Viscous  and  turbulent  terms  are  also 
adjusted  within  the  input  file.   In  summary,  the  physical  domain  is  input,  matched  to  the 
boundary  conditions,  and  transformed  to  the  computational  domain  where  the  equations  are 


solved.   The  code  was  altered  to  accommodate  the  rotating  boundary  condition,  compiled, 
and  run  on  a  Silicon  Graphics  work-station.     Specifics  of  the  code  and  input  /  output  file 
format  are  presented  in  Appendix  A. 

2 .  Grid  Generation 

A  cylindrical  31x31x41  grid  with  a  length  scale  equivalent  to  one  disk  radius  was 
developed  using  GRIDGEN .    As  flowfield  characterization  was  the  initial  purpose,  the 
complexity  of  defining  the  blade  geometries  was  not  undertaken  at  this  stage  since  it  would 
have  required  an  overlapped-grid  technique.   Therefore,  the  physical  problem  solved  was  that 
of  a  planar  rotating  disk,  with  an  axisymmetric  solution.   A  hyperbolic  grid  distribution  was 
applied  in  the  axial  direction  for  fine  resolution  of  the  disk  boundary  layer  and  swirl 
development.    Similarly,  a  hyperbolic  distribution  of  grid  lines  was  applied  radially  in  order 
to  resolve  the  flow  near  the  rim.   The  grid  domains  were  blocked  in  GRIDGEN  such  that  the 
I,J,K  directions  were  consistent  with  a  cylindrical  coordinate  system  and  an  OVERFLOW  J,K,L 
notation.   The  grid  correlation  is  given  in  Appendix  A,  Table  Al.   The  grid  is  shown  in 
Figure  1. 

3 .  Boundary  Conditions  and  Computational  Setup 

Accurately  defining  the  inlet  and  rotating  plane  boundary  conditions  was  critical. 
Although  the  code  required  an  inflow  Mach  number  to  avoid  singular  calculations,  it  was 
found  experimentally  that  Min  =  0.015  gave  a  converged  solution.   Therefore,  the  flow 
effectively  was  stationary  and  the  inflow  Mach  number  was  assumed  to  have  negligible  effect 
on  the  results.  To  effect  the  radially-varying  velocity  distribution  at  the  disk,  alteration  of  the 
code  was  required.   The  viscous,  no-slip  boundary  condition  subroutine  was  edited,  imposing 
tangential  velocities  on  the  grid  points  defining  the  disk  surface.   This  change  required  the 
upstream  flow  to  match  the  velocity  at  the  disk,  which  therefore  acted  as  a  strong  source  term 
in  the  computation.    Density  extrapolation  occurred  prior  to  the  velocity  extrapolation  and 
therefore  provided  accurate  density  updates.     Varying  the  rotational  frequency  required 
editing  the  scalar  in  the  boundary  condition  and  recompiling  the  code.    Additional  notes  on 
the  boundary  conditions  can  be  found  in  Appendix  A,  including  the  edited  subroutines.    The 
solution  method  was  fully  viscous  and  turbulence  was  estimated  using  a  Baldwin-Lomax 
turbulence  model  (Buning,1995). 

4 .  Computational  Results 

Solutions  were  obtained  for  four  test  cases  corresponding  to  the  disk  rotational  speeds 
tested  experimentally.  Table  1  summarizes  the  computational  solutions  obtained  and  lists  the 
corresponding  similarity  parameters.   Figure  2  represents  a  streamline  trace  from  the  333  Hz 


solution  and  exemplifies  the  overall  flow  structure.   Figures  3  through  6  detail  the  normalized 
static  pressure  and  normalized  stagnation  temperature  distributions  for  the  166  Hz,  266  Hz, 
285  Hz,  and  333  Hz  cases  respectively.    The  three  dimensional  figures  display  the 
axisymmetric  flowfield,  its  contours  in  the  radial  plane  and  the  solution  on  the  disk. 


Frequency 
(f) 

RPM 

(N) 

Wheel  Speed 
(U) 

Reynold's  No 
(Re) 

Mach  No 

(M) 

166  Hz 

10,000  rpm 

521.5  fps 

1.654e6 

0.442 

266  Hz 

16,000  rpm 

835.7  fps 

2.650e6 

0.709 

285  Hz 

17,000  rpm 

895.4  fps 

2.839e6 

0.755 

333  Hz 

20,000  rpm 

1046.2  fps 

3.317e6 

0.865 

Table  1.   Computational  Solutions  and  Similarity  Parameters 


The  solutions  in  Figures  3  through  6  show  the  entire  disk  and  one  peripheral  grid  location 
beyond  the  rim  to  accentuate  maxima  found  at  the  disk  rim.   Convergence  history  data  are 
given  in  Appendix  A.   Close  agreement  was  found  between  the  Mach  number  (0.442)  at  the 
disk  rim  at  166  Hz  to  that  calculated  incompressibly  using  the  method  of  Schlichting. 
Evident  in  these  solutions  is  the  small  pressure  variation  exhibited  across  the  disk  in  all  cases, 
in  contrast  to  strong  stagnation  temperature  changes.    For  the  ambient  laboratory  temperature 
of  69°F,  the  effective  recovery  temperature  on  the  disk  surface  was  estimated  to  range  from 
100°-150°F.     The  computed  pressure  ratios  of  .98  to  1.01  reflect  the  inefficiency  of  the  flat 
disk  as  a  pump  and  justify  the  use  of  blading  as  small  shock  generators  to  generate  a 
measurable  event. 


Figure  1.  Computational  Grid  (31  x  31  x  41) 
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Figure  2.  Streamline  Trace  of  Row  About  Rotating  Disk 
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Figure  3a.  —  Distribution  at  166  Hz  Disk  Frequency 
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Figure   3b.    -=—  Distribution  at  166  Hz  Disk  Frequency 
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Figure  4a.  -5—  Distribution  /  266  Hz  Disk  Frequency 


Stagnation  Temperarnrt 
1.2069 


Figure  4b.    -=—  Distribution  /  266  Hz  Disk  Frequency 
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Figure  5a.  -5—  Distribution  /  285  Hz  Disk  Frequency 
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Figure  5b.    -=—  Distribution  /  285  Hz  Disk  Frequency 
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Figure  6a.  -5—  Distribution  /  333  Hz  Disk  Frequency 
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III.     EXPERIMENTAL  SETUP  AND  CONSIDERATIONS 

Figure  7  shows  the  experimental  setup.  An  air-powered  turbine  was  used  to  rotate  the 
dual-blade  disk-rotor.   The  rotor  was  continuously  and  evenly  illuminated  and  the  steady, 
relative  pressure  distribution  was  measured  by  integrating  a  series  of  low  light  image  captures, 
with  the  final  image  representing  an  averaged  intensity  distribution.    A  once-per- re  volution 
pulse  produced  by  a  LED  /  photodiode  pair  triggered  the  image  acquisition  system  through  a 
phase-locked,  monostable  pulse  circuit.   The  acquisition  system  consisted  of  a  charge- 
coupled  device  (CCD)  camera  and  control  unit  (CCU),  operating  with  a  video  frame  grabber  / 
image  processing  system.  A  function  generator  provided  triggering  for  "wind-off  and  dark 
current  imaging.    A  description  of  the  system  components  follows. 

A.        HAMILTON  STANDARD  TURBINE-DRIVEN  FUEL  PUMP 

To  attain  a  transonic  condition  at  the  tip  of  the  disk-rotor,  a  Hamilton  Standard 
turbine-driven  fuel  pump,  Model  TPC-13,  was  mounted  within  an  existing  spin  chamber. 
Powered  by  a  300  psi,  8000  cu  ft  air  supply,  the  pump  was  capable  of  providing  wheel  speeds 
to  28,500  rpm  continuously  (Hamilton  Standard,  1968).    Turbine  mounting,  bearing 
lubrication,  air  supply,  and  disk-rotor  alignment  were  the  primary  considerations  in  the 
installation  and  are  discussed  separately. 

A  box  cradle  was  designed  and  built  from  3/4  inch  6061-T6  aluminum  plate  to 
mount  the  turbine  to  a  shelf  inside  the  protective  spin  chamber.   The  turbine  and  mount  are 
shown  in  Figures  8  and  9.   The  spin  chamber  described  by  Varner  (1995)  was  adapted  to 
accommodate  the  turbine  and  cradle  and  the  chamber  assembly  was  secured  to  a  heavy  steel 
mounting  table  inside  the  Gas  Dynamics  laboratory.   Air  supply  and  bearing  lubrication 
plumbing  were  provided  to  the  turbine  housing.   Figure  10  shows  the  housing  and  turbine  in 
place.    The  mechanical  drawings  for  the  turbine  mount  are  given  in  Appendix  B. 

Due  to  the  high  rotational  speeds,  turbine  bearing  lubrication  was  critical.   Reviewing 
lubrication  requirements  of  similar  high  performance  type  bearings  and  referring  to  the 
drawings  for  the  bearing  geometry  (Hamilton  Standard,  1968),  a  DN  value  of  1,050,000  mm- 
rpm  was  determined.   This  allowed  a  circulating  oil  lubrication  system  to  be  applied 
(Torrington  Company).    Although  the  turbine  was  designed  to  use  diverted  aviation  type  fuel 
as  the  bearing  lubricant,  military  specification  MIL-L-7808  lubricating  oil  was  substituted. 
Circulation  oil  pressure  was  provided  by  a  flushed  and  converted  hydraulic  pump,  capable  of 
providing  oil  pressures  to  3000  psi.  Because  this  far  exceeded  the  turbine  design  limitations, 
careful  control  of  the  oil  pressure  to  the  bearings  was  observed.   A  maximum  pressure  of  60 
psi  was  applied  which  forced  oil  through  the  bearing  housing  to  the  sump  return  line.   A 
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temperature-sensing  thermocouple  was  installed  in  the  return  line  prior  to  the  sump  to 
monitor  bearing  temperature.    Figure  1 1  shows  the  lubrication  pump  setup  and  Table  2 
summarizes  the  turbine  operating  limits  observed  in  the  present  testing  (Hamilton  Standard, 
1968). 


Parameter 

Operating  Limit 

Operating  Speed 

28,500  rpm1 

Outlet  Oil  Temperature 

175°F 

Inlet  Oil  Pressure 

60  psia 

Note  (1):   A  bearing  limitation.   The  disk-rotor  was  stress  limited  to  20,000  RPM. 
Table  2.   Turbine  Bearing  Lubrication  Limitations 

The  air  supply  was  plumbed  to  the  turbine  using  two-inch  steel  piping  and  a  flexible 
coupling.    A  threaded  connecting  flange  required  to  mate  the  air  supply  plumping  with  the 
turbine  air  inlet  was  manufactured  from  stainless  steel  and  secured  with  a  Marman  clamp. 
The  detailed  drawings  for  the  part  are  given  in  Appendix  B.   The  air  supply  was  controlled 
using  a  regulator  at  the  control  panel.    A  lock  nut  was  incorporated  to  constrain  the  air  supply 
valve  once  test  speed  was  attained,  thereby  limiting  speed  variation  during  the  test. 

Following  installation,  the  disk-rotor  was  aligned  to  the  turbine  drive  flange.   A 
machine  indicator  with  a  ±0.001  inch  resolution  was  clamped  to  the  mounting  bracket  for 
reference.   Axial  runouts  at  the  disk  rim  and  mid-radius  locations  were  measured  and  shim 
washers  were  placed  between  the  disk  and  drive  flange  to  bring  the  axial  rim  runout  to  within 
0.005  inches.   Table  3  lists  the  shims  required  at  each  attachment  hole,  numbered 
counterclockwise  facing  the  turbine  attachment  flange. 


f  Of  2     s   4    J  O  1 

Attachment 
Point 

Alignment 
Requirement 

1 

0.005" 

2 

0.015" 

3 

0.000" 

4 

0.003" 

Table  3.    Shim  Requirements  for  Turbine  Alignment  by  Attachment  Hole 
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B .  DISK-ROTOR 

The  test  disk-rotor,  pictured  in  Figure  12,  was  built  and  tested  previously  by  Varner, 
(1995).   It  consisted  of  a  .25  inch  thick,  6061  aluminum  disk  with  two  opposing  2.18"W  x 
.125"H  blades  designed  to  produce  a  measurable  pressure  event.   The  disk  was  adapted  to  fit 
the  turbine  drive  flange  for  the  present  study.   A  structural  analysis  was  undertaken  and 
20,000  rpm  was  found  to  give  the  radial  stress  structural  limit.   A  summary  of  the  analysis  is 
presented  in  Appendix  E.   The  radial  and  tangential  stress  distributions  based  on  the  solution 
for  a  disk  with  a  central  hole  (Baumeister,  1967),  are  shown  in  Figure  El.  A  MSC 
NASTRAN  i  modal  analysis  of  the  disk-rotor  was  also  completed.  Plate  elements  were  used  in 
the  formulation.   The  disk  hub  was  constrained  in  the  three  translational  directions  and  in  the 
transverse  and  lateral  rotational  planes.     Figure  E2  depicts  the  NASTRAN  mesh  geometry 
and  Table  El  lists  the  natural  frequencies  and  the  corresponding  modes.   The  results  show 
modes  with  complementary  frequencies  due  to  the  symmetry  of  the  model.   The  primary 
bending  and  flapping  modes  were  shown  to  occur  at  low  frequencies.   A  pitching  mode, 
however,  was  found  at  313  Hz  and  turbine  operation  in  the  vicinity  of  this  frequency  was 
monitored  for  excessive  vibration.   Figure  E3  shows  normalized  deformation  contours  for  the 
primary  modes  and  exemplifies  the  solution  symmetry  in  the  second  and  third  bending 
modes. 

C.  ONCE-PER-REVOLUTION  (1/REV)  TRIGGERING 

The  1/rev  triggering  circuit  is  shown  in  Figure  13.   A  light  emitting  diode  (LED) 
(type  SLA580LT3F)  excited  a  planar-diffused  PIN  photodiode  (PIN  3CD)  through  a  3/32 
hole  in  the  disk-rotor.  An  analog  voltage  pulse  was  generated  through  a  basic  photo- 
conductive  circuit  (UDT  Sensors,  1996).   The  pulse  was  relayed  first  to  an  operational 
amplifier  and  then  to  a  wave-shaping  circuit  to  produce  a  trigger  for  the  CCD  camera  and  to  a 
counter  to  record  the  wheel  frequency.     LED  and  photodiode  characteristics  are  provided  in 
Table  4.     Simplified  radiometry  (radiance)  and  photometry  (luminance)  conversions 
(Roberts,  1993)  were  applied  to  ensure  the  pulse  width  and  amplitude  supplied  to  the  wave- 
shaping  circuit  for  triggering  were  sufficient  over  the  operational  frequency  spectrum. 
Additionally,  the  response  of  the  photodiode  circuit  was  set  to  a  lower  order  (0.01  psec)  than 
the  gate  speed  of  the  camera  (0.1    psec). 
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Model 

Characteristic 

Specification 

LED  SLA580LT3F 

Luminous  Emission 
Wavelength  Intensity 
FWD  /  REV  Voltage 

2000  mcd 
660  nm 

2.5  V  /  4  V 

PIN3CD 
Photodiode 

Responsivity 
Capacitance  (@  28V) 
Breakdown  Voltage 
Rise  Time  (2.2RLQ) 

0.35  amps  /  Watt 
8.5  pf 
50  V 
37.4  nsec 

Table  4.    1/REV  Trigger  Component  Specifications 

The  LED,  photodiode  and  associated  wiring  were  cased  in  cylindrical  housings 
(electrically  non-conductive  for  the  photodiode  casing)  and  fit  into  a  mounting  bracket, 
shown  in  Figure  14.   The  bracket,  machined  from  .375  inch  aluminum  plate,  was  attached  to 
the  turbine  mount.   Adjustment  set  screws  in  the  bracket  allowed  axial  translation  of  the  LED 
and  photodiode  casings.   Mechanical  drawings  are  given  in  Appendix  B. 

D.        IMAGE  ACQUISITION  SYSTEM 

1 .        CCD  Camera  and  Control  Unit 

The  images  were  acquired  with  an  intensified,  electronically-gated  Xybion  ISG-350 
CCD  camera  and  its  associated  CCU  (Varner,1995).  The  camera  was  fitted  with  a  75mm,  fl.4 
Cosmicar™  television  lens  and  Oriel™  interference  filter  #53590.   The  camera  was  focused 
manually  and  a  10mm  extension  for  the  lens  was  required  for  the  present  application.   The 
camera  was  secured  to  a  tripod  which  allowed  fine  adjustment  of  the  camera  position. 
Characteristics  of  the  CCD  camera  are  given  in  the  ISG-350  users  manual  (Xybion  Electronic 
Systems,  1992).  The  interference  characteristics  of  the  filter  are  shown  in  Figure  15. 

The  8-bit  camera  dynamic  range  was  maximized  by  setting  the  CCU  camera  control 
switch  to  AUTO,  which  set  a  38mV  black  level  and  OdB  camera  gain  (Xybion  Electronic 
Systems,  1990).   This  gain  setting  avoided  potential  image  noise  created  by  increased  gain 
settings.   However,  part  of  the  8-bit  resolution  was  lost  due  to  heat-generated  noise  within  the 
CCD  elements,  which  required  the  black  level  to  be  set  at  the  noise  band  threshold.   Because 
full  video  potential  was  represented  from  0-714mV,  a  portion  of  the  video  resolution  was  lost. 
The  intensifier  gain  control  on  the  CCU-01  was  set  manually  to  a  maximum.   Due  to  short 
gate  durations,  discussed  in  the  following,  this  setting  amplified  the  low  intensities  'seen'  by 
the  camera  prior  to  their  distribution  into  the  CCD  array,  and  eliminated  the  need  for 
cumbersome  intensifier  gain  measurements.   In  addition,  the  lens  aperture  was  set  full  open  to 
the  f=1.4  setting,  maximizing  the  amount  of  light  introduced  to  the  camera  intensifier. 
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To  make  full  use  of  the  video  resolution  present,  pixel  distortion  in  the  integrated 
image  was  eliminated  by  gating  the  camera  to  freeze  the  movement  of  the  disk-rotor  blade  to 
within  one  pixel  of  travel.     Table  5  lists  the  gate  speed  requirements  for  varying  spatial 
resolutions  and  angular  frequencies  at  the  limiting  horizontal  video  resolution  of  752  pixels. 
The  gate  speed,  a  function  of  the  rotational  speed  of  the  disk,  was  set  manually  using  the 
CCU-01.  The  manual  gate  adjustment  was  made  while  displaying  the  signal  on  a  digital 
oscilloscope  via  the  gate  readout  BNC  connector  on  the  CCU-01  rear  panel. 

2 .  Illumination  System 

Paint  illumination  was  produced  using  a  1000- Watt  quartz  tungsten-halogen  lamp, 
Model  66200,  equipped  with  an  F/0.7  lens  and  controlled  by  its  associated  Oriel  lamp 
controller,  Model  6405.   The  lamp  was  filtered  using  an  ORIEL  blue-gel  filter  #66228  and 
ORIEL  interference  filter  #57521  yielding  a  380nm  excitation  source.    Lamp  controller 
voltage  was  maximized,  sacrificing  filament  lifetime  for  maximum  illumination  intensity. 
Illumination  source  and  filtering  characteristics  are  described  by  Seivwright  (1996). 
Documented  lamp  ramp-up  times  (Oriel  Corporation,  1992)  were  conservatively  observed 
prior  to  image  acquisitions. 

3 .  Frame  Grabbing  and  Image  Display 

An  EPIX  4MEG  Video  Model  12  integrated-circuit  board  and  EPIX  4MDP  V2.8 
software  were  installed  in  a  120MHz  clock  speed  Pentium™  personal  computer.   Video 
resolution  was  set  to  752H  x  480V  pixels  which  configured  the  video  memory  into  eleven 
buffered  segments.    Each  image  was  downloaded  to  the  frame-grabber  board  after  the  camera 
image  integration  period.   The  dynamic  range  of  the  frame-grabber  board  was  adjusted  to 
match  that  of  the  camera  (EPIX  Incorporated,  1993).   Acquired  images  were  displayed  for 
post-processing  on  a  Sony  PVM  1954Q  color  video  monitor.   The  monitor  allowed 
red/green/blue  (RGB)  inputs  from  standard  RS-170  video  input.   Monitor  channel  A  /  RGB 
was  used  (Sony  Corporation,  1993).   At  the  conclusion  of  the  present  work,  the  requirement 
for  increased  video  memory  was  identified  and  the  frame-grabber  board  was  upgraded  to 
EPIX  4M12-IM-64/430-14.3  (64  megabytes  of  memory)  which,  at  the  same  resolution, 
segmented  into  178  video  buffers.    Additionally,  the  software  upgrade  to  EPIX  4MIP  for 
WINDOWS95™  V3.2  was  installed. 

4  .        Monostable  Wave-Shaping  Circuit 

The  amplified  triggering  pulse  produced  by  the  LED  /  photodiode  /  op-amp  circuit 
was  converted  to  a  5V  active-low  camera  trigger  through  a  previously  developed  wave- 
shaping  circuit  (Varner,  1995).   BNC  connections  were  made  as  shown  in  Figure  7. 
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This  circuit  consisted  of  a  waveshaper,  voltage  divider  /  invertor,  TTL  logic  control  gate,  and 
monostable  pulse  circuit.    Figure  16  details  the  pulse  shaping  for  each  circuit  element  and 
Figure  17  shows  the  image  acquisition  instrumentation  setup.   The  waveshaper  and  voltage 
divider/invertor  converted  the  amplified  pulse  to  a  5V  high  TTL  pulse.      This  was  input  both 
to  the  TTL  logic  gate  and  to  a  frequency  counter  for  wheel  frequency  indication.    The  TTL 
pulse  at  the  TTL  logic  gate  was  relayed  to  the  camera  only  on  external  initiation  from  a 
Wavetek  function  generator  pulse  to  the  monostable  pulse  circuit.   The  monostable  pulse 
circuit  generated  two  5V  high  TTL  pulses,  widths  of  which  were  dependent  on  variable  RC 
time  constants.   The  first  was  the  image-inhibit  period,  during  which  the  camera  was  triggered 
and  images  integrated.    The  second  allowed  the  completed  integrated  image  to  be 
downloaded  to  the  frame  grabber.  The  monostable  pulse  circuit  was  modified  from  its 
previous  configuration  to  generate  longer  time  constants,  up  to  35  seconds.    Additionally,  a 
binary-coded  decimal  (BCD)  counter  was  built  to  count  the  trigger  pulses  relayed  to  the 
camera.   Appendix  C  contains  circuit  specifications  for  these  devices. 

5 .        Image  Capture  and  Processing 

To  best  use  the  8-bit  video  resolution  available,  image  integration  times  were  adjusted 
so  intensities  in  the  integrated  images  were  three  decimal  significant  digits  (100-255  grey  scale 
value).   This  allowed  greater  pressure  resolution  in  the  ratioed  image,  as  arithmetic  operations 
in  the  EPIX  4MIP  image  processing  software  used  the  computer  32  bit  accuracy.     It  was 
determined  experimentally  that  20.0  seconds  was  a  necessary  integration  period  for  the  gate 
speeds  required.    Image  blooming  during  this  integration  period  was  not  observed  as  reported 
by  Varner  (1995)  or  suggested  by  Xybion  (1992).    "Wind-on",  "wind-off,  dark-current, 
and  calibration  image  captures  were  scripted  with  the  available  scripting  macro  in  the  4MIP 
software  (Epix  Incorporated,  1993)  .  This  allowed  image  capture  to  be  initiated  at  the 
personal  computer  interface  and  provided  'hands-off  image  capture.    The  user  was  then 
freed  to  initiate  image  triggering  and  to  monitor  the  system.  Because  two  video  fields  were 
measured  by  the  camera,  the  default  field  value  of  0  in  the  TRIGGERED  IMAGE  CAPTURE 
menu  was  set  to  2.   Images  were  saved  using  a  file  naming  construct  given  in  Table  6.  The 
scripted  programs  for  image  capturing  and  post-processing  are  listed  in  Appendix  D. 


1st  digit 

2nd  digit 

3rd,  4th,  5th  digit 

6th,7th  digit 

suffix 

X  -  experimental 
C  -  calibration 

N  -  "wind-on" 
F  -  "wind-off 
D  -  dark  current 

###  -  frequency  (Hz) 

Test# 

.TIF 

Table  6.  File  Naming  Structure 
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Figure  7.   PSP  Experimental  Setup 
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Turbine  /  Cradle  Front  View 


Turbine  /  Cradle  Left  Side  View 


Turbine  /  Cradle  Oblique  View 


Turbine  /  Cradle  Top  View 

Figure  8.    Four- View  Aspect  of  Hamilton-Standard  TPC13  Turbine  Mount 
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Figure  10a.   Spin  Chamber  Side-View  with  Protective  Front  Plate  Installed 
(Air  Supply  Plumbing  and  1/REV  Trigger  Circuit  Wiring  shown) 


Figure  10b.  Spin  Chamber  Views; 
Turbine  Uninstalled  (Top)  and  Installed  (Bottom) 
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Temperature  Thermocouple 


Oil  Temperature  Sensor  and  Readout 


Pressurized  Supply  Li.ie  Pressure  Regulation  Valve 

Pressure  Repeater 


Figure  11.  Turbine  Lubrication  System 


24 


Disk-Rotor  Blade  Expanded  View 


Figure  12.   Dual-Bladed  Disk-Rotor 
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Figure  13.    Triggering  Circuit  Diagram 
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Wavelength 


Figure  15.   ORIEL  Camera  Filter  #53590  Characteristics 
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A.  300  mV  1/REV  Trigger 


B.  4V  Amplified  1/REV 

Trigger 


C.         9V  Active  Low  1/REV 


D.  5V  High  TTL  1/REV 


E.  5V  Trigger  Pulse 


F.  Signal  5V  Active  High 

Monostable  Pulse 


Image  Inhibit  5V  Active 
High  Monostable  Pulse 


H.  5V  Low  TTL  Signal 
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Figure  16.   Waveform  Shaping  Schematic 
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IV.  EXPERIMENTAL  PROCEDURE 
A .    ROTOR  PREPARATION 

The  disk-rotor  was  prepared  in  the  manner  described  by  Varner  (1995).    Seivwright 
describes  in  detail  the  paint  application  process  (Seivwright,  1996).   In  summary,  a  white, 
glossy  Krylon  (#1501)  undercoating  was  applied  and  allowed  to  dry  completely.   The 
pressure-sensitive  paint,  PtEOP  probe  molecules  immersed  in  a  GP-197  silicon  matrix  binder, 
was  applied  with  an  airbrush  over  the  blade  area.    A  vertical  reference  mark  was  made  at  the 
center  of  the  blade  with  a  fine  razor  for  alignment  of  "wind-on"  and  "wind-off  images. 
The  spatial  resolution  of  this  reference  mark  was  generally  2-3  horizontal  pixels  when 
imaged.  The  disk  was  attached  at  the  four  attachment  points  with  steel  bolts,  with  the  balance 
shims  (Table  2)  inserted  about  the  bolts  between  the  turbine  attachment  disk  and  the  rotor 
disk.   The  attachment  fasteners  were  torqued  to  30  in-lbs  with  a  torque  wrench,  ensuring  a 
consistent  plane  of  rotation.   The  machine  indicator,  discussed  in  the  disk  balance  section,  was 
used  to  verify  that  the  rotor  disk  rim  runout  remained  within  0.005  in.   Figure  10b  shows  the 
installed  painted  disk.  The  spin-chamber  face  plate  was  then  secured  and  the  lamp  and 
camera  set  in  place.  When  viewed,  the  blade,  at  the  twelve  o'clock  position,  was  in  the  upper 
100  vertical  pixels  of  the  display. 

B  .        CAMERA  AND  LAMP  PLACEMENT 

Upon  setting  the  lamp  and  camera  in  place,  a  spirit  level  and  inclinometer  were  used 
to  verify  the  camera  and  light  source  axes  were  normal  to  the  disk-rotor  plane  to  within  ±1.0°. 
The  camera  was  required  to  view  the  blade  at  a  4.5°±0.1°  inclination  from  the  horizontal  at  a 
distance  of  17.5  ±  0.1  inches  due  to  the  spin  chamber  face  plate  configuration.   The  lamp 
lens  assembly  was  placed  15.2  ±0.1  inches  from  the  disk  center.  The  lamp  intensity  was  set 
to  a  maximum  with  the  variable  intensity  adjustment  and  was  turned  off  between  tests  using 
the  power  switch  on  the  lamp  controller,  to  ensure  the  same  intensity  was  produced  for  both 
"wind-on"  and  "wind-off  images. 

C.        FOCUS  AND  GATING 

With  the  camera  and  lamp  positions  verified,  the  blade  edge  was  focused  manually 
with  the  camera  lens  focus  adjustment.  This  process  was  accomplished  with  the  lamp 
illuminating  the  painted  surface  in  laboratory  darkness.   The  total  paint  illumination  time  was 
monitored  due  to  the  photo-degradation  characteristics  of  the  paint.    Using  the  PIXEL  PEEK 
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/  POKE  feature  in  the  4MIP  software,  the  blade  was  measured  in  horizontal  pixels  and  its  2.18 
inch  physical  length  was  related  as  a  percentage  of  overall  horizontal  video  resolution.   The 
correct  gate  speed  from  Table  5  was  then  selected  and  set  using  a  separate  Wavetek  function 
generator  to  provide  a  simulated  trigger  (see  Figure  7). 

D.  "WIND-ON"  ACQUISITION 

The  first  measured  image  was  the  "wind-on"  acquisition.   With  the  camera  and 
circuit  settings  ready,  the  lubrication  system  for  the  turbine  was  started  and  air  was  introduced 
slowly  to  the  turbine  to  initiate  a  steady  turbine  run-up.   Oil  temperature,  atmospheric 
pressure  and  laboratory  temperature  were  recorded  during  run-up.    Once  test  operating 
frequency  was  achieved,  the  lock  nut  for  the  butterfly  valve  on  the  air  control  console  was 
secured  and  frequencies  were  maintained  to  within  ±1  Hz.  The  lamp  was  turned  on,  allowed 
to  ramp  to  its  set  intensity,  and  the  disk  was  illuminated  for  at  least  ten  seconds  prior  to 
imaging  to  allow  the  paint  time  to  respond.  "Wind-on"  script  initiation  presented  the  4MIP 
triggered-image  capture  menu  and  the  image  was  triggered  with  a  single  function  generator 
pulse.   Following  the  20.0  second  image  integration  period,  the  count  from  the  BCD  counter 
was  recorded  manually  while  the  script  cycled  to  the  next  available  buffer  and  prompted 
another  image  capture.   This  process  was  repeated  for  the  eleven  available  image  buffers, 
which  were  subsequently  averaged  with  the  4MIP  IMAGE  PROCESSING  -  IMAGE 
SEQUENCE  OPERATIONS  menu.  The  resulting  image  was  saved  in  a  .TIFF  picture  format 
in  the  hard  drive  memory  of  the  PC  under  directory  C:/4MIP/IMAGES/.   The  turbine  was  then 
shut  down  and  the  recorded  trigger  count  was  totalled. 

E.  DARK  CURRENT  ACQUISITION 

Dark  current  was  recorded  in  laboratory  darkness  with  the  camera  lens  cap  off, 
accounting  for  the  influence  of  the  LED.   A  second  Wavetek  generator  was  used  to  produce  a 
pulse  trigger  similar  to  that  produced  by  the  1/RFV  circuit  and  then  input  into  the  operational 
amplifier.    The  remainder  of  the  circuitry  settings  remained  identical  to  the  "wind-on" 
settings.   The  same  procedure  was  then  followed,  with  image  counts  recorded  over  the  eleven 
image  buffers.  The  averaged  image  was  saved  as  above.   Because  dark  current  builds  as 
circuitry  components  become  warmer,  dark  current  imaging  was  conducted  between  "wind- 
on"  and  "wind-off  to  approximate  an  averaged  'warm'  condition  of  respective  system 
components. 
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F.  "WIND-OFF"  ACQUISITION 

The  "wind-off  imaging  process  followed  the  identical  format  of  the  dark  current 
imaging.   However,  the  image  was  registered  to  the  "wind-on"  image  and  the  lamp  was 
switched  on.    This  was  accomplished  by  loading  the  "wind-on"  image  into  buffer  1  and 
marking  the  vertical  reference  mark  with  the  cross  hairs  from  within  the  PIXEL  PEEK  /  POKE 
menu.  The  CCU  gating  knob  was  set  to  AUTO  without  adjusting  the  vernier  gate  timing 
knob,  allowing  for  a  return  to  the  exact  gate  setting  used  in  the  previous  imaging  processes. 
The  lamp  was  turned  on  and  the  image  displayed  in  real  time  due  to  the  AUTO  setting 
feature.   The  disk  was  manually  moved  until  the  vertical  reference  was  located  at  the  marked 
location  of  the  "wind-on"  image.    The  gate  selection  was  reset  to  its  manual  gate  position. 
Script  initiation,  image  capture,  buffer  averaging,  and  file  saving  were  accomplished  as 
described  above. 

G .  POSTPROCESSING 

Scripting  of  the  post  processing  sequence  loaded  "wind-off,  "wind-on",  and  dark 
current  images  respectively  into  the  first  three  buffers.   The  dark  current  was  subtracted  with 
the  4MIP  IMAGE  PROCESSING  /  TWO  IMAGE  ARITHMETIC  menu  using  the  MOD  256 
subtraction  option  (Epix  Incorporated,  1993,  p. 28-1).    With  the  dark  current  subtracted,  the 
images  were  ratioed,  "wind-off  to  "wind-on",  in  the  same  menu  under  RATIO  IMAGES. 
A  ratio  coefficient,  CO,  was  selected  which  appropriately  shifted  the  ratioed  images  to  an 
optimum  resolution.    Linear  pseudocoloring  of  the  image  provided  visual  enhancement  of 
the  intensity  distribution.   The  RGB  lookup  function  is  presented  in  Figure  18.  This  process  is 
discussed  in  detail  by  Seivwright  (1996).  Pressure  ratios  were  then  estimated  through 
application  of  the  Stern- Volmer  relationship  to  the  ratioed  intensities. 
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Figure  18.     Plot  of  Linear  RGB  Pseudocoloring  Functions  and  Color  Palette 


34 


V.     RESULTS  AND  DISCUSSION 

Ratioed  images  were  acquired  for  four  rotor  frequencies;  166  Hz,  266  Hz,  285  Hz, 
and  333  Hz.     The  Io/I  ratios  were  multiplied  by  a  factor  of  50,  providing  a  common 
reference  and  allowing  continuous  greyscale  resolution.     The  technique  developed  for  this 
application  resulted  in  Table  7,  which  summarizes  each  case  and  its  image  gating,  integration 
timing,  and  average  image  counts.   The  variation  in  the  integration  periods  resulted  in  less 
than  0.63%  difference  between  "wind-off,  "wind-on",  and  dark-current  on-times.    The 
colorized  images  are  shown  in  Figures  19-22.   Image  quality  was  affected  by  image 
integration  times  during  which  noise  was  introduced  proportionally  and  resulted  in  a 
'blotchy'  effect  in  the  images.   This  effect  was  able  to  be  reduced  by  increasing  the  number 
of  images  used  in  the  averaging.  However,  the  number  of  averaged  images  was  limited  due  to 
the  photodegradation  effect  of  the  paint.   For  each  case,  the  total  PSP  excitation  time  was 
approximately  fifteen  minutes,  keeping  the  response  of  the  PSP  to  within  90%  of  its  initial 
value  (Kavandi,  1990). 

Figure  19,  the  166  Hz  case,  shows  a  nearly  constant  intensity  distribution  at  a  ratio  of 
approximately  1.2.    Proper  "wind-off  /  "wind-on"  registration  was  verified  by  the  lack  of 
blade  spatial  detail  and  constant  intensity.     Figure  20,  the  266  Hz  case,  shows  similar  traits  but 
with  the  intensity  ratio  increasing  from  1.75  to  2.25  radially.   A  vertical  reference  mark  made 
on  the  blade  can  be  identified  at  the  central,  upper  disk  rim  in  a  low  intensity.  This  case  image 
was  also  characterized  by  PSP  shedding  in  the  area  below  the  blade  trailing  edge.  This  effect 
resulted  from  horizontal  and  vertical  reference  marks  in  the  paint,  made  using  a  scalpel,  along 
which  PSP  was  shed  progressively  during  the  "wind-on"  acquisition.    Because  the 
referencing  caused  PSP  shedding,  the  referencing  was  not  used  at  higher  speeds.   The  285  Hz 
case,  Figure  21,  shows  a  similar  intensity  distribution  radially,  but  with  slightly  increased 
intensity  ratios,  from  1 .9  to  2.4.  Figure  22,  imaged  at  333  Hz,  showed  the  results  of  the 
increased  wheel  speed.   A  dramatic  intensity  change  in  the  radial  direction  was  apparent,  with 
high  intensities  shown  under  the  blade  edge.   Ratios  from  2.6  to  4.3  were  present  in  this 
image.   Also  characterizing  this  image  was  PSP  shedding  at  approximately  20%  blade  chord. 
No  reference  marking  had  been  done  in  this  region. 

Using  the  Stern- Volmer  relation  to  relate  the  intensity  maps  of  Figures  19-22  to 
surface  pressure  does  not  give  physically  meaningful  results  since  pressure  ratios  ranging 
from  1 .2  to  6.0  would  be  calculated  for  the  four  cases  based  on  previous  calibrations  of  the 
PSP  used  (Varner,  1995).   The  intensity  distributions  are,  however,  qualitative.   Figure  22, 
with  a  rim  Mach  number  of  0.865  shows  increased  intensity  under  the  blade  edge,  indicative 
of  the  stagnation  effect  on  the  radially  pumped  fluid  mass.   Because  PSP  was  shed  in  the 
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Case 

1 

2 

3 

4 

f{Hz} 

166 

266 

285 

333 

N  {rpm} 

10,000 

16,000 

17,000 

20,000 

U  {fps} 

521.5 

835.7 

895.4 

1046.2 

Mach  No 

.442 

.709 

.794 

.865 

Blade  Spatial  Resolution  (C) 

.675 

.701 

.773 

.718 

Integration  Period  (T)  {sec} 

19.87 

19.75 

20.16 

20.13 

Gate  Speed  (G)  {fisec} 

0.65 

0.40 

0.33 

0.30 

No  Images  in  Average 

11 

11 

11 

11 

Avg  Wind-Off  Images  in  Integration 

3298 

5247 

5738 

6725 

Avg  Wind-On  Images  in  Integration 

3302 

5232 

5739 

6682 

Avg  Dark  Current  Images  in  Integration 

3292 

5259 

5748 

6711 

Avg  Wind-Off  On-Time  /  Image  { (isec } 

2144 

2099 

1894 

2017 

Avg  Wind-On  On-Time  /  Image  { |isec } 

2146 

2093 

1894 

2004 

Avg  Dark  Current  On-Time  /  Image  { |_isec } 

2139 

2104 

1897 

2014 

Patm{inHg} 

29.94 

29.95 

30.03 

30.03 

Tatm{°R} 

529 

529 

529 

530 

Table  7.   Imaging  Data  Comparison  by  Case 


36 


vicinity  where  a  shock  might  begin  to  develop  in  a  transonic  flow,  the  edge  at  which  the  paint 
was  shed  may  represent  a  shock  structure.   However,  the  pressure  sensitivity  of  the  probe 
molecule  was  a  lesser  contributor  to  the  intensity  distributions  than  the  temperature  sensitivity. 
An  approximate  physical  verification  of  increased  plate  temperature  was  made  following  each 
test  by  placing  a  contact  thermocouple  against  the  disk  after  the  turbine  stopped  windmilling. 
The  increased  disk  temperatures  measured  were  consistent  with  the  numerical  stagnation 
temperature  solutions  in  Figures  3-6. 

To  examine  the  intensity  change  radially,   Figures  23a-b  graphically  compare  the 
intensities  along  two  vertical  lines  of  pixels,  in  the  radial  direction  at  15%  and  50%  chordwise 
locations.   Figure  23c  provides  a  comparison  along  a  horizontal  line  of  pixels,  the  chordwise 
direction  at  83%  radial  location.   Figures  23a-b  show  the  effect  of  radially  increasing  velocity 
on  the  local  stagnation  temperature,  which  was  verified  in  the  numerical  stagnation 
temperature  predictions.   The  increase  in  slope  between  each  case  confirms  the  dependence 
on  both  temperature  and  pressure,  as  greater  effects  are  found  at  higher  wheel  speeds. 
Figure  23c  verifies  the  uniformity  of  pressures  and  temperatures  at  nearly  constant  radius. 
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10O125LSD175    21O22S25D27S3a)3251!DJ75«O«543D475 

Vertical  Pixel  Location  {Disk  Rim  —  >  Disk  Hub} 


Figure    23a.  Vertical  IJl  Distribution  at  Horizontal  Station  X  =  15%  Blade  Chord 


75     100l25tS0     175    300225250275300325350J75.l00425450    475 

Vertical  Pixel  Location  {Disk  Rim  — >  Disk  Hub} 


Figure    23b.    Vertical  I0/I  Distribution  at  Horizontal  Station  X  =  50%  Blade  Chord 
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Horizontal  Pixel  Location   {Leading  Edge  — >  Trailing  Edge} 


Figure   23c.   Horizontal  IJi  Distribution  at  Vertical  Station  Y  =  83%  Disk  Radius 
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VI.     CONCLUSIONS  AND  RECOMMENDATIONS 

The  potential  use  of  pressure  sensitive  paint  as  a  pressure  mapping  technique  in 
transonic  turbomachinery  design  validation  has  been  furthered  by  the  successful  development 
and  application  of  a  phase-locked  imaging  technique,  providing  optimum  resolution  and 
eliminating  pixel  blur.   Incorporation  of  a  BCD  counter  allowed  verification  of  equivalent 
'on-times'  between  respective  images  while  trigger  signal  amplification  allowed  reliable 
circuit  triggering  throughout  the  frequency  range  tested.    Development  of  the  amplifier  and 
waveshaping  circuit  given  in  Appendix  C  would  enhance  the  overall  reliability  of  the  phase- 
locked  circuit,  eliminate  incompatible  circuit  components,  and  improve  system  response. 

The  turbine  installation  and  successful  operation  to  20,000  RPM  provided  a  working 
test-bed  for  further  PSP  development  in  the  transonic  range.   The  1/rev  detection  and 
associated  circuitry  proved  suitable  for  the  application.   Bearing  lubrication,  however,  was  not 
adequate  for  continuous  operation  at  higher  wheel  speeds  because  the  lubricant  was  not 
scavenged  from  the  bearing  case,  creating  high  bearing  temperatures.    Suction  pressure  on 
the  return  line  will  alleviate  this  problem.   To  allow  higher  wheel  speeds  and  the  generation  of 
a  greater  pressure  differential  on  the  disk,  a  redesign  of  the  disk  for  constant  stress  at  a  design 
condition  of  30,000  RPM  was  carried  out  and  is  presented  in  Appendix  E.   Machining  of  the 
redesigned  disk- rotor  was  initiated  at  the  conclusion  of  the  present  study. 

The  relationship  of  the  intensity  change  to  radially  changing  stagnation  temperature 
has  been  shown,  providing  the  understanding  needed  to  obtain  fully  quantitative  results. 
Because  "wind-on"  imaging  took  place  under  a  radial  temperature  change  and  "wind-off 
imaging  took  place  subsequent  to  turbine  shut  down,  the  exact  effect  of  temperature  can  not 
be  deduced.   To  effectively  quantify  the  effect  of  temperature  on  the  PSP,  the  effect  of 
oxygen  quenching  due  to  pressure  would  have  to  be  eliminated  and  only  the  effect  of 
temperature  measured.    Therefore,  using  two  "wind-on"  images,  one  with  the  effect  of 
pressure  and  one  with  the  effect  of  pressure  and  temperature,  the  images  could  in  principle  be 
reduced  using  calibration  data. 

To  obtain  quantitative  measurements  in  the  NPS  Turbopropulsion  Laboratory 
transonic  compressor  test-rig,  further  technique  development  first  must  occur  on  the  test  bed 
developed  in  the  present  study.     In  addition,  practical  PSP  application  techniques  for  the 
transonic  rotor  must  be  considered.    Recommendations  include: 

1 .  Apply  suction  pressure  to  the  turbine  lubrication  return  to  reduce  bearing 
temperatures. 

2.  Align  and  balance  the  new  disk-rotor  and  operate  to  30,000  RPM. 
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3.  Redesign  the  waveshaper  to  ensure  optimum  circuit  response  at  higher  wheel 
speeds. 

4.  Identify  an  oxygen  impermeable  binder  in  which  PtEOP  may  be  suspended 
thereby  yielding  a  temperature  sensitive  paint.   Alternatively  examine  the 
potential  of  sealing  the  current  PSP  from  oxygen,  eliminating  the  quenching 
created  by  collisions  with  oxygen.   Consideration  in  this  regard  must  be  made 
to  the  effect  of  a  sealant  on  image  distortion. 

5 .  Re-derive  the  quantum  equations  including  temperature  as  a  variable  to 
appropriately  determine  how  to  register  temperature  sensitive  and  pressure 
sensitive  images  for  quantification. 

6.  Identify  the  effect  of  blade  curvature  and  angled  camera  placement  on 
intensity  distributions. 

7 .  Determine  the  adherence  qualities  of  the  PSP  at  higher  wheel  speeds,  with  and 
without  a  Krylon  undercoating. 

8 .  Quantify  the  effect  of  the  white  Krylon  undercoating  as  an  image  enhancer 
and  determine  its  necessity  in  application  to  the  transonic  rotor,  where  paint 
thickness  and  weight  present  additional  balance  and  flowfield  distortion 
considerations. 

9.  Develop  a  registration  technique  for  "wind-off  /  "wind-on"  images  of  the 
transonic  compressor  test  rotor  accounting  for  the  mechanical  deformation  of 
the  "wind-on"  condition. 

10.  Identify  image  postprocessing  filtering  requirements  to  aid  in  reducing  signal 
noise  resulting  from  image  integration. 
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APPENDIX  A. 
OVERFLOW  NUMERICAL  SOLUTION 

A.        CODE  INPUT 

The  numerical  solutions  obtained  for  the  rotating  flat-plate  disk  were  solved  on  the 
computational  domain  previously  described,  shown  in  Figure  1.     Table  Al  describes  the 
relation  of  coordinate  systems  in  defining  the  geometry.     The  output  of  GRIDGEN  ,  the 
geometry  generating  code,  was  a  standard  PLOT3D  format.    Buning  (1995)  provides 
PLOT3D  grid  format  detail.     The  grid  was  converted  to  an  unformatted  binary  file  for  use  by 
the  OVERFLOW  solver.   The  input  grid  was  defined  as  boundary  conditions  through  the 
OVERFLOW  input  file,  'overflow.in'.   Table  A2  details  boundary  condition  definition  for 
the  Euler  and  viscous  cases.   To  verify  the  grid  compatibility,  boundary  condition 
applicability,    and  test  the  solution  independence,  an  Euler  solution  was  conducted  by 
defining  a  low  input  Mach  number  (to  M=0.21)  at  the  input  plane  and  defining  the  exit  plane 
as  an  adiabatic  wall.  This  solution  displayed  the  predicted  symmetric,  radial  flow  at  the  wall 
with  a  converged  solution  indicated  in  the  L2-norm  and  density  residuals.   A  similar 
stationary,  viscous  solution  was  computed. 


Domain                  Physical                GRIDGEN 

OVERFLOW 

Computational 

Coordinates                 6,  R,  Z                      I,  J,  K 

J,K,L 

^>C>V 

Table  Al.   Grid  Domain  and  Coordinate  Relation 


Boundary  Condition  Description 

Type 

Direction  Assigned 

Periodic 

10 

+l;J=l,allK,  allL 

Axis  (K  around) 

15 

+2;  all  J,  K=l,allL 

Characteristic  condition 

47 

-2;  all  J,  K=Kmax,  all  L 

Supersonic/Subsonic  inflow/outflow 

32 

+3;  all  J,  all  K,  L=l 

Inviscid,  adiabatic  wall  (pressure  extrapolation) 

1 

-3;  all  J,  KL2    L=Lmax 

Viscous,  adiabatic  wall  (pressure  extrapolation) 

5 

-3;allJ,K=ltoKedge3,L=Lmax 

!Euler  solns  ~  all  K 

2Viscous  solns  ~  K=Kedge  to  Kmax 

3Viscous  solns  only 

Table  A2.  OVERFLOW  Boundary  Conditions 
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Appendix  A 


Initial  solutions  were  obtained  on  the  Naval  Postgraduate  School  CRAY  processor 
through  batch  submission.   Due  to  administrative  waits  for  processor  access  and  the  relatively 
small  grid  size  used  in  the  solution,  the  OVERFLOW  code  was  compiled  on  a  Department  of 
Aeronautics  and  Astronautics  Silicon  Graphics  (SGI)  workstation.   The  application 
represented  the  first  successful  compilation  of  the  code  on  an  NPS  workstation.   Use  was 
made  of  ' makefile. sgi_mips2'  in  compiling  the  code.    The  'Makefile'  compiler  allowed 
individual  code  modification  for  deriving  the  solution  for  the  rotating  disk  at  varying  speeds. 

B .        SOLUTION  METHOD 

To  obtain  the  solution  for  a  rotating  planar  disk,  the  OVERFLOW  program  code  was 
adapted  in  the  boundary  condition  subroutine.     Table  A3  gives  the  path  and  subroutine  calls 
which  were  modified,  in   descending  hierarchy. 


.../overl  .6aw/ns/bc.F 

.../overl  .6aw/ns/bc/bcvwpx.F 

.../overl  .6aw/ns/bc/bcnosl.F 

Table  A3.  Modified  OVERFLOW  Subroutines 

The  rotation  of  the  disk  was  produced  by  imparting  a  radially  varying  velocity  distribution  in 
V  and  'w'  at  each  grid  point  which  defined  the  viscous  disk.   In  producing  this  no-slip 
rotating  frame,  the  actual  code  modification  required  was  minimal.   The  primary 
modification  was  made  in  bcnosl.F.   First  the  physical  grid  coordinates  were  entered  as 
arguments  in  the  respective  subroutine  call.   The  coordinates  were  required  to  apply  the 
appropriate  radial  location  for  tangential  velocity  computation  from  the  rate  of  rotation. 
Second,  the  internal  call  to  subroutine  uvwrho.F,  which  provided  time  metrics  for  the 
contravariant  velocities,  was  commented  out.  Next,  the  appropriate  velocity  distribution  was 
computed  and  subsequently  multiplied  by  the  density.   Because  updates  to  the  density  (qj) 
occurred  prior  to  the  velocity  extrapolation  at  the  next  time  step,   correct  updates  were  made 
to  the  solution  vector  Q,  shown  by  Equation  Al. 

[q,  q2  q3  q4  q5]  =  [p  pu  pv  pw  e]  <A1> 

To  obtain  the  appropriate  velocity  distribution  at  the  disk  surface,  the  internal  call  to 
bczero.F  was  retained  to  impart  a  zero  V  (axial)  velocity.   The  V  and  'w'  velocities  were 
then  computed  from  the  tangential  velocity,  Ve,  in  non-dimensional  form.     Additional  inputs 
to  the  solution  file  included  the  implicit  time  step,  'dT',  a  minimum  Courant-Freidrichs-Levy 
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(CFL)  number,  and  turbulence  parameters.       Buning  (1995)  provides  recommendations  for 
these  solution  variables.   The  parameters  listed  in  the  'overflow. in'  file  provided  remained 
constant  for  all  solution  cases. 

The  calculation  of  the  non-dimensionalized  V  and  'w'  velocity  components  at  the 
disk  boundary  follows  in  Equations  A2-A3.   Non-dimensional  quantities  are  shown  with 
asterisks.  Taking  wheel  speed  in  RPM  (Q),  the  angular  rotational  rate  was  given  by 


^  =  Qtc  (A2). 


30 


The  angular  frequency  was  non-dimensionalized  by  the  freestream  speed  of  sound,  a^,  and 
the  length  scale  of  the  grid;  one  disk  radius,  R. 


co*=  -^i  (A3)" 


The  velocity  components  in  the  standard  Cartesian  sense  were  transformed  from  the  non- 
dimensional  tangential  velocity,  Ve*,  which  was  obtained  from  multiplication  of  the  angular 

frequency,  go*,  by  the  radial  location,  f .   The  axisymmetric-to-Cartesian  transformation  was 
made  by  incrementing  the  angle  by  quantity  d0  and  resolving  the  v*  and  w*  components 
based  on  the  angular,  or  'J',  position.    The  'JKLmax'  subscripts  generalize  the  location  on 
the  disk  for  which  the  computation  was  solved.  Equations  A4-A7  show  the  velocity 
computation. 

V*  =  coY  (A4). 


(A5). 

(A6). 
(A7). 


de  = 

J 

2tc 
-1 

max 

JKL 

Ve' 

'sin(J* 

d9) 

max 

w* 

JKL 
max 

Ve' 

'cos(J*d9) 
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C.  SOLUTION  OUTPUT 

The  solutions  represented  in  Figures  3  through  6  were  graphically  portrayed  using 
the  computational  fluid  dynamics  solution  modeler  FAST  .   The  solution  set  was  output  in 
unformatted  binary  construct  in  a  file  'q.save'.   This  file  was  read  into  FAST  for  graphical 
display.     Additional  output  files  included  the  residual  data  in  'resid.out'.   This  file,  in  text 
format,  was  used  to  observe  the  convergence  history  of  the  solution  sets.  Figure  Al  shows  the 
L2-norm  and  density  residual  history  for  the  333  Hz  case.   This  convergence  history  is 
representative  of  all  the  cases  presented  in  this  study.   2500  iterations  proved  sufficient  in 
attaining  converged  solutions. 

D.  PROGRAM  CODES 

The  following  codes  are  listed  for  reference.   The  'Overflow. in'  file  was  used  in  each 
test  case  for  consistency.   Modifications  made  within  the  OVERFLOW  subroutines  'bcnosl.F', 
'bc.F',  'bcvwpx.F'  are  denoted  in  boldface  type  and  are  remarked  in  italics. 

1 .  Overflow.in 

2.  bcnosl.F 

3.  bc.F 

4.  bcvwpx.F 
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L-2  Norm  and  Density  Residual  Convergence  History 


1000  1250  1500 

Numerical  Iteration 


2000 


2500 


Figure  Al.  333  Hz  Case  Convergence  History 
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1.   'OVERFLOW.IN' 

$GLOBAL 

CHIMRA=  .F.,    NSTEPS=2500,    RESTRT=  .F.,    NSAVE  =2500, 

NQT    =0, 

SEND 
SFLOINP 

ALPHA  =0.0,  FSMACH=  0.015,  REY    =  6.00E6,  TINF  =520.000, 

SEND 
SVARGAM 

SEND 
SGRDNAM 

NAME  =  'Viscous  Flow  Against  Disk,  31x31x41', 
SEND 
SNTTERS 

SEND 
SMETPRM 

IRHS  =0,    ILHS  =2,    IDISS  =2, 

SEND 
STIMACU 

DT     =  0.5,  ITIME=  1,  TFOSO  =  1.00,  CFLMIN=0.01, 

SEND 
SSMOACU 

ISPECJ=2,  DIS2J  =  2.00,  DIS4J  =  0.2, 

ISPECK=2,  DIS2K  =  2.00,  DIS4K  =  0.2, 

ISPECL=2,  DIS2L  =  2.00,  DIS4L  =  0.2, 

SMOO  =  1.00, 

EPSE  =0.35, 

SEND 
SVISINP 

VISCJ  =  .F.,  VISCK  =  .F.,  VISCL  =  .F., 

NTURB=    1, 


ITTYP=    1, 

ITDIR=    -3, 

JTLS  =    1, 

JTLE  =  31, 

KTLS  =    1, 

KTLE  =    31, 

LTLS  =    1, 

LTLE  =  41, 

TLPAR1=  0.3, 

SEND 

SBCINP 

NBC    =    6, 

IBTYP  =     5, 

1, 

47, 

47, 

14, 

10, 

IBDIR=    -3, 

-3, 

3, 

-2, 

2, 

1, 

JBCS  =        1, 

1, 

1, 

1, 

1, 

1, 

JBCE  =       31, 

31, 

31, 

31, 

31, 

1, 

KBCS  =       1, 

21, 

1, 

31, 

1, 

1, 

KBCE  =      20, 

31, 

31, 

31, 

1, 

31, 

LBCS  =      41, 

41, 

1, 

1, 

1, 

1, 

LBCE  =      41, 

41, 

1, 

41, 

41, 

41, 

SEND 

SSCEINP 

SEND 
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APPENDIX   B. 
TURBINE  APPARATUS  MECHANICAL  DRAWINGS 

1 .  Hamilton  Standard  TPC13  Turbine  Mounting  Bracket  Sideview  /  Topview 

2.  Hamilton  Standard  TPC13  Turbine  Mounting  Bracket  Front  /  Back  View 

3.  Hamilton  Standard  TPC13  Turbine  Mount  /  Spin  Chamber  Interface 

4 .  LED  /  Photodiode  Mount  Bracket  Views 

5.  Turbine  Connecting  Flange 
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Figure  B4.   Light  Emitting  Diode  /  Photodiode  Bracket  Detailed  Drawing 
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Figure  B5.   Turbine  Connecting  Flange  Detailed  Drawing 
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APPENDIX   C. 
CIRCUIT  DESIGN 


1 .  BCD  Counter  Circuit 

2 .  Power  Supply  Circuit 

3 .  Waveshaper  Redesign  Proposal 
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APPENDIX   D. 
4MIP  TOOL  SCRIPTED  MACROS 

The  scripting  of  the  imaging  process  allowed  for  repetitive,  'hands  off  measurements 
and  gave  the  tester  the  ability  to  monitor  the  overall  process.   The  scripts  were  developed 
using  the  TEACH  /  EXECUTE  SCRIPT  menu  in  the  EPIX  4MIP  TOOL  software.  This  menu 
option  recorded  keyboard  stroke  inputs  in  a  executable  text  file.    Subsequent  execution  of 
the  text  file  allowed  the  software  to  emulate  the  recorded  keystrokes.   Because  of  the  large 
memory  requirement  per  image,  the  software  was  best  able  to  acquire,  store,  and  manipulate 
the  data  by  first  executing  the  SPLIT  SCREEN  option  to  alleviate  the  need  to  exit  and  switch 
between  menus  during  the  script  execution  (EPIX  Users  Manual,  1993).   The  scripts  were 
saved  under  the  path  and  suffix  \4MIP\SCRIPTS\name.MPX. 

The  imaging  process  commenced  with  the  execution  of  the  script  PSP_PIC.MPX 
which  first  executed  a  split  screen  and  commenced  a  series  of  image  captures  to  fill  the 
available  1 1  segmented  memory  buffers.   A  sub-script  denoted  CHECKER.MPX  was  then 
initiated  three  times  in  a  sequence.   This  script  execution  resulted  in  a  video-like  play  back  of 
the  1 1  buffers,  providing  a  visual  check  for  error  within  the  data.    With  the  data  checked,  the 
script  then  executed  another  sub-script,  PSP_AVG,  which  averaged  the  1 1  images  in  the 
SEQUENCE  OPERATIONS  menu.   The  resulting  image  was  then  saved  as  IMAGE1.TIF. 
The  process  was  then  repeated  with  the  second  image  saved  as  IMAGE2.TIF.  The  final  image 
was  developed  by  loading  and  averaging  the  two  images.  The  file  was  saved  using  the  file 
naming  protocol  shown  in  Table  6.   The  scripts  listed  in  this  Appendix  follow: 

1 .  PSP_PIC.MPX 

2.  CHECKER.MPX 

3.  PSP  AVG.MPX 
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APPENDIX  E. 
DISK  STRESS  ANALYSIS  AND  ROTOR  REDESIGN 

The  structural  evaluation  of  the  flat  plate  disk-rotor  consisted  of  stress  and  modal 
analyses.  The  structural  analysis  determined  the  limiting  rotational  speed  of  the  disk  and 
verified  characteristic  modes  of  the  disk  would  not  be  excited  at  the  wheel  speeds  used  in  the 
experiment.    The  stress  analysis  was  based  on  strain-displacement  relationships  for  a  uniform, 
thin  rotating  disk.   Aerodynamic  stresses  and  rotor-blade  stresses  were  not  incorporated  in  the 
computation.   The  modal  analysis  used  finite-element  methods  in  the  structural  analysis  code 
MSC-NASTRAN.   A  summary  of  each  analysis  is  provided.   Due  to  the  limiting  speed  of  the 
current  disk- rotor,  the  maximum  tip  Mach  number  attainable  was  0.865.   To  obtain  sonic  tip 
conditions,  a  new  disk-rotor  design  is  required.   A  proposed  design  of  this  disk  for  use  on  the 
PSP  test-rig  apparatus  is  presented  and  summarized. 

A.        STRESS   ANALYSIS 

From  the  development  of  the  strain-displacement  relationships  for  a  uniform,  thin 
rotating  disk  with  a  central  hole  ,  the  following  relationships  for  radial  and  tangential  stress 
was  obtained  (Baumeister,  1967): 


3+1)       2/r>2 ,  1-D 

par(R2+ 

4    K  3+\) 


a9  =  i-pipa)2(R2+4_V  )  (El) 


Or  =  2±^pco2(R2-r2)  (E2) 

o 


For  type  6061-T6  aluminum  plate  (density  =  0.1  lbsm/in3,  yield  stress  =  40,000  psi),  the 
maximum  wheel  speed  was  limited  by  tangential  stress  at  the  inner  hole  location  to  21,750 
RPM.  The  tangential  stress  at  the  inner  radius  is  plotted  versus  wheel  speed  in  Figure  El 
(top).   Because  blade  and  aerodynamic  stresses  were  not  computed,  the  operational  limit  was 
reduced  to  20,000  RPM.   Figure  E2  (bottom)  plots  the  radial  and  tangential  stresses  at  the 
20,000  RPM  operational  limit,  as  a  function  of  the  radial  distance  from  the  rim. 
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B  .        MODAL  ANALYSIS 

The  modal  analysis  of  the  flat-plate  disk-rotor  was  undertaken  using  MSC-NASTRAN 
for  WINDOWS  finite  element  methods.    Plate  elements  were  used  in  the  formulation.  The 
disk  hub  was  constrained  in  the  three  translational  directions  and  in  the  transverse  and  lateral 
rotational  planes.    Figure  E2  depicts  the  NASTRAN  mesh  geometry  and  Table  El  lists  the 
frequencies  and  the  corresponding  modes.   The  results  show  modes  with  complementary 
frequencies  due  to  the  symmetry  of  the  model.   The  primary  bending  and  flapping  modes 
were  shown  to  occur  at  low  frequencies.  A  pitching  mode,  however,  was  found  at  313  Hz  and 
turbine  operation  in  the  vicinity  of  this  frequency  was  monitored  for  excessive  vibration. 
Figure  E3  shows  normalized  deformation  contours  for  the  primary  modes  and  exemplifies 
the  solution  symmetry  in  the  second  and  third  bending  modes. 


0)(Hz) 

Natural  Mode 

Cu(Hz) 

Natural  Mode 

CD  (Hz) 

Natural  Mode 

17.850 
17.972 

1st  Bending  1 

119.47 

2nd  Rotating 

269.90 
271.68 

7th  Bending 

21.232 

1st  Flapping2 

128.15 
129.70 

Mixed  mode5 
Bending  / 
Flapping 

30297 
303.99 

5  th  Flapping 

26.101 
27.043 

2nd  Bending 

147.90 
147.91 

5th  Bending 

312.24 
313.15 

1st  Pitching4 

56.973 

3rd  Bending 

161.94 
167.84 

3rd  Flapping 

330.60 
340.01 

Mixed  mode 
Pitching  /  Bending 

87.319 

1st  Rotating3 

205.10 
207.31 

6th  Bending 

341.93 
342.79 
343.39 

Mixed  mode 
Bending  / 
Flapping 

98.318 
98.489 

4th  Bending 

227.31 
228.31 

4th  Flapping 

IBending  modes  are  asymmetric  about  the  plane  of  rotation  (XY  plane). 

2Flapping  modes  are  symmetric  about  the  plane  of  rotation  (XY  plane). 

3Rotating  mode  is  symmetric  about  the  axis  of  rotation  (Z  axis). 

4Pitching  mode  is  translation  in  the  plane  of  rotation  (XY  plane). 

5Mixed  modes  have  characteristics  of  more  than  one  of  the  previously  defined  modes. 

Table  El.  Disk-Rotor  Modal  Frequencies  (Hz) 
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C.        DISK-ROTOR  REDESIGN 

The  design  for  a  disk-rotor,  capable  of  bearing  the  stress  associated  with  sonic  tip 
speeds,  is  based  on  constant  stress  distribution.  The  constant  stress  design  requires  varying 
the  disk  thickness,  which  subsequently  will  increase  the  difficulty  in  the  disk  manufacturing. 
Applying  radial  equilibrium  and  assuming  constant  stress,  equation  E3  is  derived  (Mattingly, 
1987). 


3.  =  e   2°v        R 


(E3) 


To  obtain  an  estimate  of  disk  thickness  at  the  material  yield  and  tensile  stress  limits,  Figure  E4 
is  plotted.   The  physical  requirement  for  mounting  the  disk-rotor  to  the  turbine  is  a  flat 
interface,  however,  resulting  in  the  estimated  design  geometry  portrayed.    In  addition, 
imaging  requirements  at  the  disk  rim  necessitate  a  portion  of  the  outer  rim  be  flat.   The 
design  drawing  is  presented  in  Figure  E5. 
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Figure  El.  Disk-Rotor  Stress  Distributions 
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V:  Disk  Rotor    ZX  Front  View 
C:  Constrained 


Y 


V:  Disk  Rotor  XY  SidaView 
C:  Constrained 


Y 

L. 


Figure  E2.  NASTRAN  Disk-Rotor  Finite  Element  Model 
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1st  Bending  Mode 


1st  Flapping  Mode 


3rd  Bending  Mode 
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2nd  Flapping  Mode  1st  Pitching  Mode 

Figure  E3.  Disk-Rotor  Modal  Characteristics 
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